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Abstract 

We consider a relativistic strongly interacting Bose gas. The interaction is 
manifested in the off-shellness of the equilibrium distribution. The equation of 
state that we obtain for such a gas has the properties of a realistic equation of 
state of strongly interacting matter, i.e., at low temperature it agrees with the 
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one suggested by Shuryak for hadronic matter, while at high temperature it 
represents the equation of state of an ideal ultrarelativistic Stefan-Boltzmann 
gas, implying a phase transition to an effectively weakly interacting phase. 
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1 Introduction 

Quantum chromodynamics (QCD), the fundamental theory of strong interactions, 
can be perturbatively solved only in the region of asymptotic freedom, i.e., for large 
momenta of quarks and gluons [|TJ. At low momenta, quarks and gluons interact 
strongly and are confined inside hadrons. In this case, the expansion parameter of 
perturbation theory, the running coupling constant a s , is of the order of unity, so 
that perturbative methods are not applicable and one has to use non-perturbative 
methods. One such method is the construction of effective theories to describe the 
behaviour of hadronic matter, such as the a — u model . However, one of the most 
successful non-perturbative methods is the use of lattice gauge calculations ||, which 
are especially suitable for studying perturbative as well as non-perturbative effects 
in QCD. The presently available lattice data mostly concern simulations of SU(N) 
pure gauge theory |4], ||, since the treatment of dynamical fermions on a lattice is 
difficult ||. Moreover, lattice artefacts are believed to be well under control only in 
the pure gauge case ||. A striking feature of lattice simulations of SU(2) and SU(3) 
pure gauge theory is a phase transition (of apparently first order for SU (3) [£|, |], [?j 
and second order for 577(2) theory ||) from a phase of confined gluons ( "glueballs" ) 
to one of deconfined gluons ("gluon plasma"), leading to a sharp rise in the energy 
density of a gluon g function of temperature at a phase transition temperature 

T c @j H (Figs. 1,2). In the a — uj model, a similar phase transition (or, at least, a 
rapid increase of the energy density in a small temperature interval) appears at small 
net baryon densities || (Fig. 3). This is due to a strong enhancement of the scalar 
meson interaction at T c ~ 200 MeV, leading to a transition from a phase of massive 
baryons to a phase of massless baryons. Apparently, a phase transition to a weakly 
interacting phase seems to be a fundamental feature of strongly interacting matteif] 







To understand the lattice data from a simple physical point of view, Rischke et 
al. 1 10, have constructed a phenomenological model for the gluon plasma, in 



x It should be, however, mentioned that a phase transition in the a — uj model may be an artefact 
of the approximations made in the calculation of the thermodynamic quantities in Fig. 3 (in the 
mean- field approximation) pl|. 
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which gluons with large momenta are considered as an ideal gas with perturbative 
corrections of order 0(a s ), while gluons with low momenta are subject to confining 
interactions and do not contribute to the energy spectrum of free gluons. The equation 
of state for this model, although it quantitatively reproduces the lattice data for 
the thermodynamic functions of £77(3) pure gauge theory above the deconfinement 
transition temperature (Fig. 2), it has a rather complicated form and is therefore not 
suitable for practical use, e.g., in astrophysics (for description of stellar structure), 
or in cosmolody (for treatment of a hadron-plasma phase transition in the early 
universe). In this paper, we suggest another equation of state, which reflects the 
main properties of strongly interacting matter (i.e., a phase transition to a weakly 
interacting phase) and agrees qualitatively with the lattice data for £77 (2) pure gauge 
theory both above and below the transition temperature, and has a much simpler form, 



as compared to that of ref. [10 



First, let us note that in strongly interacting matter, particles undergoing con- 
tinual mutual interaction are necessarily off-shell. Therefore, the effect of strong 
interaction in such a system may be represented by the off-shellness of its particles. 
The equilibrium state of such a system should be characterized by a well-defined rel- 
ativists mass distribution around the on-shell value. Thus, instead of dealing with 
interaction explicitly, we reduce the problem to the description of the relativistic 
off-shell ensemble. The role of interaction then consists in determining the effective 
thermodynamic parameters governing the mass distribution in a strongly interacting 
system. 

The physical framework for the description of a relativistic off-shell ensemble has 
been established by Horwitz and Piron [|13] as a manifestly covariant relativistic dy- 



namics whose consistent formulation is based on the ideas of Fock |L4j and Stueck- 
elberg |0|, in which the four components of energy- momentum are considered as 
independent degrees of freedom, permitting fluctuations from the mass shell. In 
this framework, the dynamical evolution of a system of N particles, for the classical 
case, is governed by equations of motion that are of the form of Hamilton equa- 
tions for the motion of N events which generate the space-time trajectories (particle 
world lines) as functions of a continuous Poincare-invariant parameter r ]13| [15|j . 
These events are characterized by their positions g M = (t, q) and energy-momenta 

= (E, p) in an 8iV-diHiensional phase-space. For the quantum case, the system 
is characterized by the wave function Vv(<?i, (hi ■ ■ ■ > In) £ L 2 (R 4N ), with the measure 
d A qid A q2 • • • d A q N = d AN q, = gf ; /i = 0, 1, 2, 3; i = 1, 2, . . . , N), describing the dis- 
tribution of events, which evolves with a generalized Schrodinger equation [fl3[| . The 
collection of events (called "concatenation" |]16|) along each world line corresponds 
to a particle, and hence, the evolution of the state of the iV-event system describes, 
a posteriori, the history in space and time of an TV-particle system. 

For a system of N interacting events (and hence, particles) one takes |U| (we 
use the system of units in which h = c = kg = I; we also use the metric g^ v = 
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(-,+,+,+)) 

^Ew + %fe-,?N), (i.i) 

i 

where M is a given fixed parameter (an intrinsic property of the particles), with the 
dimension of mass, taken to be the same for all the particles of the system. The 
Hamilton equations are 

dq- = dK_ = pf 

dr dp^ M' 

M = _^ = _^L. (i.2) 

dr dq itJ , dq ifi ' 
In the quantum theory, the generalized Schrodinger equation 

d 

i-Q-^riqi, ?2, • • • , ?tv) = KiPAqi, q*, ■ ■ ■ , qN) (1.3) 

describes the evolution of the iV-body wave function ip T (qi, q2, ■ ■ • , (Zjv)- 

In the present paper we restrict ourselves to a relativistic Bose gas, in order to 
compare the results with experimental data of pure gauge theory lattice simulations. 
We show that our results agree with those for 577(2) pure gauge theory. It should 
be, however, noted that since the underlying theory is basically different from QCD, 
a comparison with the 577(2) lattice data can only be qualitative. From this point of 
view, the similarity is remarkable. 



2 Ideal relativistic Bose gas 

Gibbs ensembles in a manifestly covariant relativistic classical and quantum mechan- 



ics were derived by Horwitz, Schieve and Piron [17]. To describe an ideal gas of 
events obeying Bose-Einstein statistics in the grand canonical ensemble, we use the 
expression for the number of events found in |l7j (for our present purposes we assume 
no degeneracy), 

N = V«) E n*. = £ 1 (2.1) 

fcM fc M e^-wxzsi)™ - 1 

where is the system's four-volume and m 2 = —k 2 = —k^k^ is the variable 
dynamical mass. Here, in addition to the usual chemical potential //, there is the 
mass potential \i K corresponding to the Lorentz scalar function K(p,q) (Eq. (1.1)), 
here taken in the ideal gas limit, on the iV-event relativistic phase space; in order to 
simplify subsequent considerations, we shall take it to be a fixed parameter (which 
determines an upper bound of the mass distribution in the ensemble we are studying, 
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as we shall see below). To ensure a positive-definite value for n^, the number of 
bosons with four- momentum k^, we require that 



2 

m 



m-fi-fiK—^O. (2.2) 

The discriminant of Eq. (2.2) must be nonnegative, which gives 

M , . 

" £ w (2 - 3) 

For such fi, (2.2) has the solution 





For small [ifix/M, the region (2.4) may be approximated by 

fi<m<— . (2.5) 

One sees that fix plays a fundamental role in determining an upper bound of the 
mass spectrum, in addition to the usual lower bound m > /i. 

Replacing the sum over (2.1) by an integral, one obtains for the density of 
events per unit space-time volume n = N/V^ ||18|| , 



1 rm 2 m 3 dm sinh 2 3 d3 , 
n =-r^ 2 , 2.6 

where m x and m 2 are defined in Eq. (2.4), and we have used the parametrization [[TjJ 

p° = m cosh /3, 

p 1 = m sinh 3 sin 9 cos 0, 

p 2 = m sinh /3 sin 6 1 sin (f>, 

p 3 = m sinh/? cos 9, 

< 9 < 7T, < < 27T, -oo < /3 < oo. 

In what follows we shall take fi ~ (as for the case of the ensemble of gauge 
bosons). The integral (2.6) is calculated in ref. (in the high-temperature Boltz- 
mann approximation for the integrand): 



T 4 r _ 2 „ , ,i 2M 2m 



471-3 



2 - ^XoU 



x = — ee — (2.7) 



where K v (z) is the Bessel function of the third kind (imaginary argument), m c = 
Mj fix is the central value around which the mass of the particles are distributed, in 
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view of (2.4), and 2m c is the upper limit of the mass distribution in our case of small 
fi, in view of (2.5). One calculates the characteristic averages |20j to be 

8 — x 3 K 3 (x) 

{E} = r »-4) ' (2 - 8) 

2 — x 2 K 2 (x) 
,2v T 2 16-x 4 ir 4 (x)+4x 3 K 3 (x) 

" J 2-x 2 if 2 (x) ' 1 j 

(P 2 ) = 3^ = 3T<£>, (2,1) 

and obtains the following thermodynamic functions (the particle number density, 
pressure and energy density) 0: 

= (,») = II izf^W, (2 , 2 ) 



1,™, _ r* 8-x 3 - K ,( , 
\ J /yu — 2 

,oo\ T 4 40 -z 4 Ar 4 (z) +x 3 K 3 (a ) 



V = =;(T u )gu = — 2 = iV T, (2.13) 



= (T ) = - 4V (2.14) 

7r z ar 



where (J M ) and (T 7 ^) are the average particle four-current and energy-momentum 
tensor, respectively, given by fl9fl : 

(J") = ^^), ( T ^> = T -ff^V v )- (2-15) 

In (2.15), Tav is the average passage interval in r for the events which pass through 
the small (minimal typical) four-volume AV in the neighborhood of the _R 4 -point; it 
is related to a width of the mass distribution around the central value, Am, as follows 

T AV Am = 2tt. (2.16) 

The expressions (2. 13), (2. 14) for p and p are obtained from Eq. (2.15) for (T^). 
They are, moreover, thermodynamically consistent. One may varify easily that the 
relation 

P = T %~P ( 2 - 17 ) 

is satisfied [EDI. 



2 If we had not used the Boltzmann limit for the integrand in (2.6), one would obtain the factors 
C(3) w 1.202 in Eq. (2.12) and ((4) = tt 4 /90 » 1.082 in Eqs. (2. 13), (2. 14). 
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For low T, Eqs. (2. 13), (2. 14) reduce, through the asymptotic formula 

4u 2 - 1 



1 + 



, z»l, 



to 



p 



2T 6 



7r 2 m 2 



p 



10T 6 



n 2 m 2 



5p, 



(2.18) 



(2.19) 



consistent with the "realistic" equation of state suggested by Shuryak for strongly 
interacting hadronic matter P2| . 



For high T, we use another asymptotic formula |24 

~2 



KJz) 



1 - 



4(i/ 



+ 



2 << 1, 



and obtain 



P 





(2.20) 

(2.21) 
(2.22) 



where psb = 7" 4 / tt 2 and psb = 3psB are the pressure and energy density of an ideal 
ultrarelativistic Stefan-Boltzmann gas. Therefore, as T —>■ oo, the thermodynamic 
functions of the relativistic Bose gas we are considering become asymptotically those 
of a Stefan-Boltzmann gas, implying a phase transition to an effectively weakly in- 
teracting phase, from the phase of strong interactions described by Eq. (2.19). 
It follows from Eq. (2.10), via the asymptotic formulas (2. 18), (2. 20), that 



(m 2 ) 

J^2 



T « 2m r 



2m 2 c /T 2 , T » 2m c 



(2.23) 



The dependence of p/psb, p/Psb and (m 2 )/T 2 on temperature are shown in Figs. 4,5. 
At T ~ 0.2 m c (corresponding to z ~ 0.1 in Figs. 4,5), there is a smooth phase 
transition to a weakly interacting phase described by Eqs. (2. 21), (2. 22). 



3 Concluding remarks 

The manifestly covariant framework discussed in the present paper can be an effective 
tool in dealing with realistic physical systems. The equation of state (2.12)-(2.14) 
obtained in our work reflects the main properties of strongly interacting matter (i.e., 
a phase transition to a weakly interacting phase), and agrees qualitatively with the 
lattice data for SU(2) pure gauge theory. 
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The question naturally arises of why the SU(2) lattice data || appear to contain 
a second order phase transition but the SU(3) data [|| |5| appear to contain first order 
one. There is no a priori difference between SU(2) and 577(3) pure gauge theories 
from a statistical mechanical point of view. As shown in ref. ||, SU(3) pure gauge 
theory simulations on 24 3 x lattices indicate a rapid rise in p + p as a function 
of temperature which takes place in the case of = 4, reflecting a sharp first order 
phase transition. This rise is broadened considerably in the case of N? = 6, with the 
slope of the curve diminished by almost a factor 3, indicating a smoother transition in 
this case. As remarked by the authors, this softening of the structure of the transition 
as Nt is increased may well continue as the continuum limit is aproached. Thus, in ref. 
Q the apparent first order nature of the transition in the case of SU(3) pure gauge 
theory has been called in question. Moreover, there are indications from lattice QCD 
calculations that when fermions are included, the phase transition may be of second 
or higher order |Z5| . Altogether, these observations suggest that a realistic equation of 
state of strongly interacting hadronic matter should be expected to contain a second 
or higher order phase transition, as reflected by the equation of state obtained in our 
work. 

As remarked by Ornik and Weiner f26fl , a single equation of state which, at high 
temperature, describes a quark-gluon phase and, at low temperature, a hadronic 
phase and which contains a phase transition of either first or higher order provides 
a more satisfactory theoretical description than one in which each phase is described 
by a different equation of state. In this way, the equation of state that we obtained 
can be considered as a candidate for such a realistic equation of state of strongly 
interacting matter, in contrast to the equations of state of, e.g., Rischke et al. [[Hj 
and Shuryak |22j each of which describes just one of the phases (above and below 
the transition, respectively). The introduction of the quark degrees of freedom in 
this equation of state, as well as taking into account an effective interaction potential 
in an explicit form in the strongly interacting phase, and perturbative corrections in 
the weakly interacting phase, should enable one to derive an equation of state which 
we expect to be more accurate for the description of the phenomena taking place in 
strongly interacting hadronic matter. The derivation of such an equation of state and 
its possible implications in astrophysics and cosmology are now being worked out by 
the authors. 
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FIGURE CAPTIONS 



Fig. 1. The energy density of SU(2) pure gauge theory on a 10 3 x 3 lattice as a 
function of T/Al. Taken from ref. [8]. 



Fig. 2. Fit of the thermodynamic functions, according to the equation of state 
of ref. [10], to the lattice data for SU(3) pure gauge theory. Taken from ref. [10]. 



Fig. 3. Phase transition in the a — to model. Taken from ref. [10]. 



Fig. 4. The pressure and energy density, according to the equation of state (2. 13), (2. 14), 
as functions of z = T/2m c . 



Fig. 5. The average mass squared, according to Eq. (2.10), as a function of 
z = T/2m c . 
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